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Abstract. We present a survey of recent results concerning a remarkable 
class of unitary matrices, the CMV matrices. We are particularly interested 
in the role they play in the theory of random matrices and integrable systems. 
Throughout the paper we also emphasize the analogies and connections to 
Jacobi matrices. 



1. Introduction 

From many points of view, Jacobi matrices play a central role among all Her- 
mitian matrices. We wish to discuss a class of unitary matrices recently introduced 
by Cantero, Moral, and Velazquez, [CanMorVellj . which play a similar role among 
unitary matrices. While this analogy may be extended in many directions, we in- 
vestigate here their occurrence in the theory of random matrices and integrable 
systems. 

By Jacobi matrix we mean an n x n tri-diagonal matrix 

"61 ai 



a„_i 
a„_i bn 

with aj > 0, bj e M. 

Given coefficients ao,...,Q;„_2 G the unit disk in C, and a„_i € S^, let 
pk — — and define 2x2 matrices 

\_Pk -oik_ 

for < fc < n — 2, while S_i — [1] and S„_i — [q;„_i] are 1x1 matrices. From 
these, form the n x n block-diagonal matrices 

£ = diag(E:o, 22,24, .. .) and X = diag(S_i, Si, S3, . . . ) . 

The CMV matrix associated to the coefficients uq, . . . ,a„_i is C = CA4. It is a 
5-diagonal unitary matrix, given by: 
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Following |Siml| . we will refer to the numbers ak as Verblunsky coefficients. 
A related system of matrices was discovered independently by Tao and Thiele, 
[TaoThij . in connection with the non-linear Fourier transform. Their matrices are 
bi-infinite and correspond to setting odd-indexed Verblunsky coefficients to zero. 
(This has the effect of doubling the spectrum; cf. jSimll p. 84].) 

Both Jacobi and CMV matrices arose in the study of orthogonal polynomials, so 
it is natural that we begin there. We will first describe the relation of orthogonal 
polynomials to Jacobi matrices and then explain the connection to CMV matrices. 

Given a probability measure dv supported on a finite subset of R, say of car- 
dinality n, we can apply the Gram-Schmidt procedure to {1 ^ tij J tij ^ m m m J tj^ ^} and 

so obtain an orthonormal basis for L^{dv) consisting of polynomials, {pj{x) : j = 
0, . . . , rt — 1}, with positive leading coefficient. In this basis, the linear transforma- 
tion f{x) ^ xf{x) is represented by a Jacobi matrix. An equivalent statement is 
that the orthonormal polynomials obey a three-term recurrence: 

xpj{x) = ajPj+i{x) + hjPj{x) + aj-ipj-i{x) 

where a„i = and p„ = 0. 

We have just shown how measures on R lead to Jacobi matrices; in fact, there 
is a one-to-one correspondence between them. Given a Jacobi matrix, J, let dv be 
the spectral measure associated to J and the vector ei — [1, 0, . . . , 0]"^. Then J 
represents f{x) t-^ xf{x) in the basis of orthonormal polynomials associated to dv. 

To explain the origin of CMV matrices, it is necessary to delve a little into the 
theory of orthogonal polynomials on the unit circle. For a more complete description 
of what follows, the reader should turn to |Siml| . 

Given a finitely-supported probability measure d^ on 5^, the unit circle in C, we 
can construct an orthonormal system of polynomials, (pk, by applying the Gram- 
Schmidt procedure to These obey a recurrence relation; however, to 
simplify the formulae, we will present the relation for the monic orthogonal poly- 
nomials $fc(z) 

(1.3) $fe+i(z)=z$fc(z)-afe$^(z). 
Here $^ denotes the reversed polynomial: 

k k 

(1.4) <^u{z)^Y.^iz' ^ $^(z) = 5] CfcV = . (i)' 

i=0 (=0 

and ak are recurrence coefficients. They are in fact the same as the Verblunsky 
coefficients which appeared above, in the definition of the CMV matrix. When 
djj, is supported at exactly n points, 6 D for < < n — 2, while Q!„_i is a 
unimodular complex number. 

The Verblunsky coefficients completely describe the measure djj.: There is a 1- 
to-1 correspondence between probability measures on the unit circle supported at 
n points and Verblunsky coefficients (ao, . . . , a„_i) with G D for < fc < n — 2 
and a„_i S S^. 

From the discussion of Jacobi matrices, it would be natural to consider the matrix 
representation of f{z) ^ zf{z) in L^{dfi) with respect to the basis of orthonormal 
polynomials. This is not a CMV matrix; rather it is a Hessenberg matrix which 
Simon, |Siml| . has dubbed a GGT matrix, from the initials of Geronimus, Gragg, 
and Teplyaev. Perhaps the most striking difference from a CMV (or Jacobi) matrix 
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is that a GGT matrix is very far from sparse — generically, all entries above and 
including the sub-diagonal are non-zero. 

Cantero, Moral, and Velazquez had the simple and ingenious idea of applying 
the Gram-Schmidt procedure to {1, z, z~^, z^, z~^, . . . } rather than {1, z, . . .}. The 
resulting functions, Xk{z) (0 < fc < n — 1), form a basis and are easily expressible 
in terms of the orthonormal polynomials. In this basis, the map f{z) i— » zf{z) is 
represented in an especially simple form: 

Theorem 1 (Cantero, Moral, Velazquez, 2003). In the orthonormal basis {xk{z)} 
of L^{d^), the operator f(z) ^ zf{z) is represented by the CMV matrix C associated 
to the Verblunsky coefficients of the measure dfi. 

The C = CM factorization presented above originates as follows: Let us write 
Xfc,0<fc<n— 1, for the orthonormal basis constructed by applying the Gram- 
Schmidt procedure to {1, z^^, z, z^^, z^, . . . }. Then the matrix elements of C and 
A4 are given by 

Cj+i.k+i = {X]{z)\zxk{z)), Mj+i.k+i = (xj(z)lxfe(z)). 

See |Siml| for further discussion. 

The measure d/i can be reconstructed from C in a manner analogous to the 
Jacobi case. Let dfi be the spectral measure associated to a CMV matrix, C, and 
the vector ei. Then C is the CMV matrix associated to the measure dfj,. 

Proofs of these statements can be found in .CanMorV cU or S imT] . 

An interesting special case is when the measure d/z is symmetric with respect 
to complex conjugation. This occurs if and only if all Verblunsky coefficients are 
real. It is a famous observation of Szego (see |Szel §11.5]) that the polynomials 
orthogonal with respect to this measure are intimately related to the polynomials 
orthogonal with respect to the measure dv on [—2, 2] defined by 

(1.5) / f{z + z-^)dii{z)^ f f{x)dv{x). 

Moreover, the recurrence coefficients for these measures are related by the Geron- 
imus relations: 

^fc+i = (1 - a2fc-i)a2fc - (1 + a2fc-i)a2fc-2 

(1.6) 1/2 

Ofc+l = {(1 - a2i;-i)(l - a2fe)(l + Q:2fc+l)} 

The classical proof of (|1.6(l relies on the relation between the orthogonal polynomials 
associated to the two measures. An alternate proof was given by Rowan Killip and 
the author, |KilNenll Proposition B.3], who derive H1.6|l from an explicit relation 
between the CMV and Jacobi matrices of the two measures dfi and dv, respectively. 

2. Circular /3-ensembles 

In 1962, Dyson |Dys| introduced three ensembles of random unitary matrices 
with a view to simplifying the study of energy level behavior in complex quantum 
systems. Earlier work in this direction, pioneered by Wigner, focused on ensem- 
bles of Hermitian matrices. In both of these situations, the eigenvalue distributions 
coincide with the Gibbs distribution for particles of Coulomb gas at inverse temper- 
ature P = 1, 2, 4 on the circle or the real line, respectively. Dumitriu and Edelman, 
[DumEde^ . constructed tri-diagonal matrix models for two of the three standard 
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examples of the Coulomb gas on the real line for any inverse temperature: the /3- 
Hermite and the /3-Laguerre. Rowan Killip and the author obtained the analogous 
result on the circle, using CMV matrices. 

One of the main features of all of these models is that they are defined in terms of 
0{n) independent random variables which appear in matrix entries via very simple 
combinations. Moreover, the matrices are all sparse, having 0(n) non-zero entries. 
Finally, these are not asymptotic models (see ^ForRail, and |Lip| ), but exact ones. 
(Throughout this section, n denotes the number of particles.) 

Let us elaborate: The easiest way to obtain a normalizable Gibbs measure on 
the real line is to add an external harmonic potential V{x) = . This gives rise 
to the probability measure 

(2.1) E(/)0C f{xi,...,Xn)\^{xu...,Xn)f\{e-^'^''^Uxi---dXn 

3 

on M", where, as usual, A denotes the Vandermonde determinant. This is known 
as the Hermite ensemble, because of its intimate connection to the orthogonal 
polynomials of the same name, and when (3 — I, 2, or 4, arises as the eigenvalue 
distribution in the classical Gaussian ensembles of random matrix theory. Dumitriu 
and Edelman showed that H2.1|l is the distribution of eigenvalues for a symmetric tri- 
diagonal matrix with independent entries (modulo symmetry). The diagonal entries 
have standard Gaussian distribution and the lower diagonal entries are 2^^/^ times 
a x-distributed random variable with the number of degrees of freedom equal to P 
times the number of the row. 

The second example treated by Dumitriu and Edelman is the Laguerre ensemble. 
In statistical circles, this is known as the Wishart ensemble, special cases of which 
arise in the empirical determination of the covariance matrix of a multivariate 
Gaussian distribution. For this ensemble, one needs to modify the distribution given 
in (|2.1|l in two ways: each particle Xj is confined to lie in [0, oo) and is subject to the 
external potential V{x) = —a log(a;)-|-a;, where a > —1 is a parameter. In DuniEde] 
it is shown that if i? is a certain nxn matrix with independent x-distributed entries 
on the main and sub-diagonal (the number of degrees of freedom depends on a, /3, 
and the element in question) and zeros everywhere else, then the eigenvalues of 
L = BB^ follow this distribution. 

The third standard ensemble on the real line is the Jacobi ensemble. The dis- 
tribution is as in (|2.1|) . but now the particles are confined to lie within [—2, 2] and 
are subject to the external potential V{x) = — alog(2 — x) ~ 61og(2 + x), where 
a,b > —1 are parameters. This corresponds to the probability measure on [—2, 2]" 
that is proportional to 

(2.2) I A(xi, ...,Xn)fl[i2- x,yi2 + x,f dxi--- dx„. 

i 

Dumitriu and Edelman did not give a matrix model for this ensemble, listing it as 
an open problem. 

On the unit circle, one does not need a confining potential in order to define the 
Gibbs distribution for n particles of the Coulomb gas at inverse temperature /3; it 
is given by 

(2.3) KU)^ j-j /(e^«^...,e^^")|A(e^^^...,e^^")|^|i...|^ 
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for any symmetric function /. CMV matrices realize a random matrix model with 
this eigenvalue distribution. 

We say that a complex random variable, X, with values in the unit disk, D, is 
6i/-distributed (for ;/ > 1) if 

(2.4) mix)} - ^ // fiz){l - Vft-^-^l^ Sz. 



For > 2 an integer, this has the following geometric interpretation: If v is chosen 
from the unit sphere 5"^ in R"^"*"^ at random according to the usual surface measure, 
then v\ + iv2 is 0i,-distributed. As a continuation of this geometric picture, we 
shall say that X is 0i -distributed if it is uniformly distributed on the unit circle in 
C. With this definition, we have the following family of matrix models: 

Theorem 2 (Killip, Nenciu, 2004). Given (3 > 0, let ak ^ 0/3(ri-fc-i)+i be indepen- 
dent random variables for < k < n — 1. Then the CMV matrix C ~ CM defined 
by these Verblunsky coefficients gives a (sparse) matrix model for the Coulomb gas 
at inverse temperature (3. That is, its eigenvalues are distributed according to the 
Cibbs distribution (|2.3|l . 

Finally, we present a tri-diagonal matrix model for the Jacobi ensemble H2.2|l . 
The independent parameters follow a beta distribution: A real- valued random vari- 
able X is said to be beta-distributed with parameters s, i > 0, which we denote by 
X - B{s,t), if 

(2.5) E{/(X)} = ^'7(,!)?ftf j\ ^^""^^^ - '^^'"'(^ + '^^^ 

Note that B{!^,^) is the distribution of the first component of a random vector 
from the i/-spherc. 

Theorem 3 (Killip, Nenciu, 2004). Given P > 0, let ak, < k < 2n - 2, be 
distributed as follows 



(2.6) ak 



Q^ 2n-k-2 p j^a^l^ 2n^2^/3 + 6+1) k CVCn, 

g^ 2n-k-3 p + a + b + 2, ^''~]^-^ l3) k odd. 



Let a2n-i — 0L-\ = —1 and define 

(2.7) 6fc+i = (1 - a2k~\)a2k - (1 + a2k~\)a2k-2 

(2.8) Ofc+i = {(l-a2fc-i)(l-a2fc)(l + a2fc+i)}'^' 
for < fc < n — 1; then the eigenvalues of the tri-diagonal matrix 

bi ai 



J 



ai 



are distributed according to the the Jacobi ensemble H2.2|l . 

Let us note that one of the essential ideas behind the proof of this theorem is 
to construct the distributions on Verblunsky coefhcients and use the Geronimus 
relations H1.6|l to transfer them to the Jacobi setting. 
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3. The Ablowitz-Ladik system: Structure 

By analogy with the connection between the Toda lattice and Jacobi matrices, 
we investigate the existence and properties of an integrable system related to or- 
thogonal polynomials on the unit circle and CMV matrices. The main evolution of 
the system is defocusing Ablowitz-Ladik (also known as the integrable discrete non- 
linear Schrodinger equation). An excellent elementary introduction into integrable 
systems is given in jDei| : [BabBerTalj is a detailed presentation of many aspects of 
the subject. For the more advanced notions that we shall use (especially regarding 
Lie and Poisson-Lie structures) we refer to |01sPerReySem| . 

The celebrated Toda lattice is the central example of a discrete integrable Hamil- 
tonian system. It models a 1-dimensional chain of particles with exponential nearest 
neighbor interactions and it was introduced by Morikazu Toda, |Tod| , to explain the 
findings of the famous 1955 Los Alamos computer experiment of Fermi, Pasta, and 
Ulam, |FerPasUla| . Complete integrability of the system was proved by Flaschka in 
1974, |Fla| ■ by introducing a change of variables that allowed him to set the system 
in Lax pair form. Indeed, if Flaschka's variables and bk are set to be the entries 
of a Jacobi matrix, J, then the evolution of the system is given by the equation 

i = [j,p], 

with 

ai 
-ai 02 



P= J+-J- = 



-02 



Thus the Toda lattice is intrinsically connected to Jacobi matrices and orthogonal 
polynomials on the real line R. 

In |Sim2[ Section 11.11], Barry Simon and the author found that the integrable 
system naturally related to the theory of orthogonal polynomials on the unit circle is 
the defocusing Ablowitz-Ladik system. It was introduced in 1975-76 by Ablowitz 
and Ladik, [AblLadll IAblLad2] . as a space-discretization of the cubic nonlinear 
Schrodinger equation (NLS). It reads 

(3.1) -«/3„ = (/3n+i + /3„-i) - 2(3,,, 

where /S^I/SnjcIDisa sequence of complex numbers inside the unit disk and 
The analogy with the continuous NLS becomes transparent if we rewrite (jH.lll as 

-iPn = Pn+1 - Wn + Pn-l ~ |/3„P(/3„+l + /^n-l)- 

From this point onwards, / denotes the time derivative of the function /. 
If we change variables to ak{t) — e^**/3fc(i), this system becomes 

(3.2) -iak = pi{ak+i + "fc-i), 

which is a little simpler. If we then choose a_i and an-i to lie on the unit circle, 
then they do not move and we obtain a finite system of ODEs for a^, < fc < n — 2, 
which is the specific case we treat. 

Nenciu and Simon, |Sim2l Sectionll.il], found that if periodic Verblunsky co- 
efficients are allowed to evolve according to Hamiltonians in the Ablowitz-Ladik 
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hierarchy, then their movement takes place on tori which are isospectral with re- 
spect to the CMV matrix. The symplectic form they introduced extends naturally 
to the case of finite Verblunsky coefficients, which we present here. In this setting 
it is given by: 



(3.3) {/,5}al 



(3.4) 

for any smooth functions /, g 
and, as usual. 



n-2 
n-2 



df dg df dg 



duj dvj 



dvj duj 



df dg df dg 



daj daj 



where a,- 



daj daj 



- ivj for all < j < n — 2, 



d 


1 / d 


. d 


da 


~ 2\du 


dv 



d 



d 



and — = -!—+ 1 — 

ou ov 



da 2 

For clarity, we will call this bracket the Ablowitz-Ladik bracket. Note that this 
differs by a factor of two from that used in ISim2| and |Nen) : we will adjust the 
results accordingly. 

On the set of Verblunsky coefficients (aoi ■ • • j ctn-2) £ B"^^ (with fixed Q!„-i G 
S^) we consider the Hamiltonians Re(iirm) and Im^Km), where 

for TO > 1. These form the Ablowitz-Ladik hierarchy. Indeed, the evolution of the 
Verblunsky coefficients under the flow generated by Jie{Ki) is the Ablowitz-Ladik 
evolution: 

{aj,Re(ifi)}AL = oij = + aj+i). 

The evolution equations for the CMV matrices under these flows in the Ablowitz- 
Ladik bracket were found by Ncnciu INenj : 

Theorem 4 (Nenciu, 2005). The Lax pairs for the to"^ Hamiltonian of the finite 
defocusing Ablowitz-Ladik system are given by 

(3.5) {C, Re{K^)}AL = [C, iC™ + i(CJ')*] 

and 



(3.6) 



{C,lmiK,^)}AL = [C,CT-iC^r] 



as the matrix with entries 





ifj < k; 




if 3 = k; 


0, 


if j > k. 



for all m > 1, where for a matrix A we set A 
iA+)jk -- 



The system of nonlinear differential-difference equations 

(3.7) q;„ = (1 - a^)(Q:„+i - q;„_i), {a^} C (-1, 1), 

is known as the discrete modified KdV equation (see [AblLadl] and [Gekp or the 
equation of the Schur flows (see |FayGek| ). This system's main interest lies in 
its connection to the Toda and Volterra (or Kac-van Moerbeke) lattices. More 
precisely, the Schur flows appear in the work of Ammar and Gragg |Amm(Traj as 
an evolution equation on Verblunsky coefficients obtained by transferring the Toda 
equation via the Geronimus relations from the a's and 6's. As such, it is an evolution 
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on real a's. One can then relax this condition and think of it as an evolution on 
complex coefficients, which preserves reality (that is, if the a's are real at time 0, 
then they remain real for all times). 

As it turns out, this description applies to "half of the Ablowitz-Ladik flows, the 
ones generated by lm(_R'm) for all n > 1. Indeed, the evolution on the Verblunsky 
coefficients under these Hamiltonians is equivalent to the Lax pairs 

{C,lmiKra)}AL = [C,(C'")+ - ((C™) + )1, 

which preserve reality of the a's. So it is natural to ask ourselves what are the 
corresponding evolutions of the a's and b's. Moreover, we recover (j3.7|l as the 
evolution generated by —lm{Ki) under the Ablowitz-Ladik bracket (j^.^fl . 1)^1 4|l . 
In fact, in the finite case, all the evolutions induced by Im(isr„) become, via the 
Geronimus relations, simple combinations of the evolutions in the Toda hierarchy 
(this can be deduced from the evolution of the corresponding spectral measures 
from |KilNen20 . 

There exists an abstract way of defining Poisson brackets on Lie algebras. In 
particular, if q denotes the (associative) algebra oi n x n complex matrices, then 
the algebra structure gives rise to a natural Lie algebra structure: 

[B, C]^BC~ CB. 

Of course, this also results from viewing q as the Lie algebra of GL{n;C). As a 
vector space, g = [ © a, where 

a = {A: A = ^A''}, 

which is the Lie algebra of the group l]{n) of n x n unitary matrices, and 

1 = {A e g : Lij = for i > j and Lis £ 

which is the Lie algebra of the group L(n) of n x n lower triangular matrices with 
positive diagonal entries. We will write tto and tti for the projections into these 
summands. 

This vector-space splitting of g permits us to give it a second Lie algebra struc- 
ture. First we define the classical R-matrix i? : g — > g by i? = ttj — tTo. The second 
Lie bracket can then be written as either 

[X,Y],^^[R{X),Y] + ^[X,R{Y)] VX,yeg,or 

[L + A,L' + A']^ = [L, L'] - [A, A'] \f L,L' e [, and A, A' e o. 

Using this R-matrix, one can define several Poisson brackets on the Lie algebra 
g. The Lie-Poisson (or Kirillov) bracket associated to the Lie algebra (g, [,]^) is 
linear, and it is relevant to Jacobi matrices: under the embedding J >—>■ iJ, the 
manifolds of Jacobi matrices with fixed trace are symplectic leaves. In our case, we 
focus on the Gelfand-Dikij (or quadratic) bracket, denoted here by {-j-Igd- We 
refer to |01sPerReySem| for the definitions and proofs of the general results of the 
theory. 

Theorem 5 (Killip, Nenciu, 2005). The manifold of CMV matrices with fixed 
determinant form a symplectic leaf in the Poisson manifold (g, {^-jcn)- 

The proof of this theorem relies on the fact a imitary matrix having the same 
shape as a CMV matrix is exactly a CMV matrix. Moreover, one can show that 
the set of CMV matrices with fixed determinant is an orbit under the action of 
a certain group of dressing transformations (see Theorem 4.7 of |KilNen2] ) . In a 



CMV MATRICES: A SURVEY 



9 



concurrent paper, Luen-Chau Li has independently derived this result. In a 
sense, his approach is the reverse of ours: he begins by investigating the action 
of certain dressing transformations, while we arrive at their existence only after 
studying the problem by other means. 

As is natural in the theory of Poisson-Lie algebras, we focus on Hamiltonians 
given by (t){B) = Imtr/(_B), where / is a polynomial. The general theory guar- 
antees that the evolution equations they generate will be expressible as Lax pairs. 
In our case these Lax pairs at a CMV matrix are the same as those obtained from 
Theorem 0] by taking the appropriate linear combinations. Given this result, the 
following is not surprising: 

Theorem 6 (Killip, Nenciu, 2005). For any < k,l < n - 2 the Gelfand-Dikij 
brackets of the Verblunsky coefficients are given by 

(3.9) {ak,ai}GD = and {ak,ai}GD = -2iSkipl. 

That is, the restriction of the Gelfand-Dikij bracket to CMV matrices agrees 
with the Ablowitz-Ladik bracket when written in the coordinates given by the 
Verblunsky coefficients. From now on we denote both brackets by {•,•}. 

4. The Ablowitz-Ladik system: Asymptotics 

Shortly after the discovery of the Lax pair representation for the Toda lattice, 
Moser, |Mosj . gave a complete solution for the finite system. Specifically, he discov- 
ered the angle variables associated to the actions of Henon 'Hen' and Flaschka |Flaj . 
In addition, he studied the long-time asymptotics of the system and determined the 
scattering map. 

The identification of the Ablowitz-Ladik bracket as an abstract Poisson-Lie 
bracket allows one to follow the same path in the unitary setting. More precisely, 
we begin by investigating the evolution under a Hamiltonian of the Ablowitz-Ladik 
hierarchy of the spectral measure M = associated to a CMV matrix C and 

the vector ei = [1,0,..., 0]-^: 

Theorem 7 (Killip, Nenciu, 2005). Under the flow generated by 4>{B) = Imtr{/(i3)}, 



(4.1) dt (log[/i,(t)]) - {0,log[M,]} = F(e^«0 - E^(e'^'')w(0, 

1=1 

where F{z) ~ 2Rcz/'(z). Consequently, 

(A 0) „ (f) - exp[F(e'^Ot]/.,(0) 

and for any j, I G {1, . . . , n — 1}, 

(4.3) {0l,^\og[^iJ/^ln]}^SJl■ 

A special case of this result (for Schur flows) has already appeared, |MukNak| . 
The approach used there was to begin with a special case of (|4.2() and determine 
the induced evolution on the Verblunsky coefficients. 

In defining fik and Zk, we can choose any ordering we please; however, there is a 
particular condition on this choice that simplifies the formulae below. Namely, we 
require that 

(4.4) Ai > A2 > • ■ • > A„, 
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where we use the shorthand Xk — F{zk) = 2 Re[zfc/'(zfc)]. Of course genericahy, the 
ordering will be strict. Note that by the continuity of F this labeling is well-defined 
on an open set which is invariant under the flow of 0. 

In what follows we will assume that we are indeed in a situation where the 
ordering described above is strict. This simplifies the formulae, but is not necessary; 
the interested reader can find the full results in KilNen2l . 

Under the flow generated by the Hamiltonian the masses have the following 
asymptotics 



(4.5) 



log[^fc(t)] = -{Xi- Xk)t + \og 



Ml 



+ 0(6-"*) 



ast —^ oo. (The quantities for which we do not specify time dependence are assumed 
to be evaluated to t = 0.) In particular ^i(i) — *■ 1 and, if fc > 1, then //^(i) 
exponentially fast. We now turn to the asymptotics of the Verblunsky coefficients: 



Fix 1 < fc < 



1. As all A," are distinct, 



(4.6) 



where 



(4.7) 



ak-iit) = (-1)'^ 'zi-'-Zk 



1 + Cfe-ie 



-(Afc — Afc-|-i)i 



0(e 



-7M 



6- 



(zfcZfc+l - 1) 



Mfe+l 
Mfe 



fe-1 

n 

1=1 



Zk+l 



Zl 



Zk - Zl 



and 7 > (Afe - A^+i) > 0. 

In other words, if all Xj are distinct, then viewed as a curve in the disk, ak~i{t) 
approaches the boundary in a fixed non-tangential direction. This simply amounts 
to the statement that ^fe-i is non-zero and arg(^fe_i) = arg(zfcZfc+i — 1) belongs to 
(— 7r/2, 7t/2). Let us also note that the asymptotics of pk-i are easily deduced from 
(Oil : 



pl_,{t) - -2Re(a-i)e-(''=-''=+^)*-fO(e-^*) 



\Zk+l — Zk 



2 Mfc+l 

Mfe 



fc-i 

n 



Zk+l - Zl 



Zk - Zl 



This shows that the factors C{t) and A^(i) of the CMV matrix C{t) diagonalize 
as t ^ 00 and hence so does C{t). Moreover, the eigenvalues are ordered by the 
corresponding value of _F(2). This is a well known phenomenon for the Toda lattice. 

Let us note here that when the Xj are not all distinct, C{t) converges to a direct 
sum of CMV matrices and their adjoints. Specifically, if 



A 



fe-1 



> Xk 



Xk 



> Xk4 



then afe_i, . . . , ak+m-2 do not approach the boundary and C(oo) has a non-trivial 
block of size m beginning at row/column k. If k is odd, this will be a CMV matrix; 
if k is even, it will the adjoint of a CMV matrix. While this phenomenon cannot 
occur for the Toda lattice, it can occur for Hamiltonians in the same hierarchy. 
Some examples of non-diagonalization are discussed on page 389 of | DeiLiTom| . 

Finally, our understanding of the asymptotic behaviour of the relevant quantities 
allows us to find canonical coordinates for the system. We know from the Lax pairs 
that the eigenvalues Poisson commute: {9j,6k} — 0, for all 1 < j, fc < n, where 9j is 
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the argument of the eigenvalue zj . Moreover, the evolution of the spectral measure 
given in Theorem 13 implies that 

{ei,^\0g[flj/Hn]} ^Sjl. 

It would be natural to imagine that 

Oi,... ,6'„_i, ilog[/ii//x„], . . . , ilog[^„_i//x„] 

are canonical coordinates. This is the case for the Toda lattice (as can be concluded 
from DeiLiNanTomJ). The same does not happen here; instead, the following holds: 

Theorem 8 (Killip, Nenciu, 2005). For any labelling of the eigenvalues, 

(4.8) {log[/i2/Mi],log[Ai3M]} =2cot(^) +2cot(^) +2cot(^) 

in the Gelfand-Dikij (or Ahlowitz-Ladik) bracket. 

While the previous results concerning asymptotics for the Ablowitz-Ladik system 
are proved using the abstract (Gelfand-Dikij) form of the bracket, we were unable 
to do the same here. Instead note that {log[/x^//Lt/], log[/Ltfe/^/]} is constant under 
the flow generated by a Hamiltonian <f> defined as above, and so one can work 
in a regime where the masses are already exponentially ordered. We then use 
perturbation theory to obtain more precise asymptotics for \og[iij / ni] for j — 2,2, 
in this regime, and plug these into the explicit (Ablowitz-Ladik) bracket to obtain 
the result. It would be interesting to find a proof using only the Gelfand-Dikij 
bracket, and to better understand the appearance of the cotangent function in 
formula (|4.8|l . 

In an earlier paper, [KilNenl] , Rowan Killip and the author considered the map 



(4.9) {Ol,fJ,l, ■ ■ ■ ,9n-l,fJ,n-l,dn) 1"^ {uo,Vo, . . . , U„_i , i;„_2 , 0) 

where at — Uk + ivk and q;„_i = e**. We found that: 

Proposition 9. The Jacobian of the change of variables 1)4. 9|l is given by 



det 



9(6*1, /^i, . . . ,6'„) 



= -2i 



/il • • • ^„ 



The analogous result for Jacobi matrices is due to Dumitriu and Edelman, 
[DumEde] . which served as our guide. Both of these proofs are probabilistic in 
nature. 

We also asked for a simpler, more direct derivation of the Jacobian. A recent 
preprint of Forrester and Rains, jForRai2| . gives a very direct, computational so- 
lution. Percy Deift showed us a derivation of the Jacobi matrix result using the 
symplectic structure naturally associated to the Toda lattice. The same can be 
done in the circular case (see jKilNen2l Corollary 8.6]). The key idea is the fol- 
lowing: As we can write the underlying symplectic form in either set of variables, 
we can view (|4.9f) as a symplectomorphism between two concrete symplectic man- 
ifolds. In particular, it must preserve the Liouville volume, which allows one to 
prove Proposition |5| 
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